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Abstract 

Dynamics of a charged relativistic particle in a uniform magnetic field and an obliquely propa- 
gating electrostatic shock wave is considered. The system is reduced to a two degrees of freedom 
Hamiltonian system with slow and fast variables. In this system, the phenomenon of capture into 
resonance can take place. Under certain condition, a captured phase point stays captured forever. 
This corresponds to unlimited surfing acceleration of the particle. The preprint is a more detailed 
version of a comment on the paper by D.Ucer and V.D.Shapiro ffl, intended for the Comments 
section of Physical Reviews Letters. 
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In Letter 0, unlimited surfing acceleration of relativistic particles by a shock wave normal 
to a uniform magnetic field was considered. The mechanisms of surfing acceleration were 
studied in many papers, mainly in the case of acceleration by a harmonic wave. The aim 
of this comment is to point out that the methods developed in [[J allow for more detailed 
study and further progress in the topic of shock wave acceleration. In particular, the case of 
an oblique shock wave can be considered. Also, it can be shown that particles with initial 
velocities far from the shock wave velocity can also be captured in the mode of unlimited 
acceleration. 

Like in fl[], consider a charged relativistic particle of charge e and rest mass m in a uniform 
magnetic field B and an electrostatic shock wave of potential $ = — $ tanh(kq — ut), where 
$o > 0, uj > 0, uj/\k\ = u is the phase velocity of the shock wave, q is the radius vector. 
Choose an orthogonal coordinate system (q\, q 2 , g 3 ) such that B = B e 3 is along the g 3 -axis 
and k lies in the (qi, g 3 )-plane, k = (hi, 0, k 3 ). The Hamiltonian function of the particle is 

H = {m 2 c A + c 2 pl + c 2 p\ + (cV 2 - eB oqi ) 2 ) 1/2 - e$ tanhfagi + A; 3 g 3 - ut), (1) 

where "P 2 = P2 + eB qi/c and p = (pi,P2,P3) is the particle's momentum. Introduce nota- 
tions: 

OJc = — , k = (kl + kl) 1/2 , e = Q c = u c /e, sina = k 3 /k. 

mc mc z 

Consider the problem in the following range of parameters: |p|/(mc) ~ 1, uj/(kc) ~ 1, e <C 

1, uj c jijj ~ e. Rescale the variables: p 13 = p 13 /(mc), gi j3 = eqi^/c, ki t s = fci )3 c, H = 

H/(mc 2 ). Following 0, one can canonically transform ([l|) into the form (tildes are omitted): 

H = -ul+ [l + k 2 (I + pcosa/k) 2 +p 2 sm 2 a + n 2 c q 2 ] 1/2 -etanh(f) = no-etanh(f), (2) 

where canonically conjugated pairs of variables are (p, e~ l q) and (J, 0) , = k±qi + k 3 q 3 — ut. 
The corresponding Hamiltonian equations of motion imply that while 0^0, total change 
in variable / is a value of order e, and hence the trajectory of the particle in the (p, q, /)- 
space lies in a vicinity of the intersection of a second order surface 7io =const and a plane 
/ =const. This intersection is an ellipse corresponding to the Larmor motion. However, 
along a trajectory that crosses the resonance <fi = &Ho/dI = the value of I can change 
significantly. 

The resonant condition dH /dI = defines a surface / = I res (p, q) in the (p, q, J)-space, 
called the resonant surface. The condition implies that projection of the particle's velocity 
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FIG. 1: Phase portraits of the "pendulum-like" system, a) b = 0.5, b) b = 1.2. 



onto the direction of vector k equals the phase velocity of the wave. Intersection of the 
resonant surface and the surface Tio = const is a second order curve whose kind depends 
on the parameter values (see [0). This curve is called the resonant curve. The motion 
in a neighborhood (of the width of order y/e) of the resonant surface possesses certain 
universal properties (H, 01)- I n particular, the Hamiltonian F = TC/e of the particle in this 
neighborhood can be written in the form: 

F = e-'Aip, q) + F (P, 0, p, q) + O(Vi), (3) 

where A(j>, q) is TCo restricted onto the resonant surface, P = (I — I res (p, q))/V^ + 0{\/e) = 
0(1), and canonically conjugated pairs of variables are (P,<j>) and (p, e~ 3 ' 2 q). The function 
F is so-called "pendulum-like" Hamiltonian, and in the case under consideration it is F = 
g(j), q)P 2 /2 — tanh</> + b(p, q)<p, where 



b(p, q) 



fl 2 . cos a 



(k 2 -tu 2 y/ 2 (l+p 2 sin 2 a + fi2 g 2)i/2' 



g(p, q) 



k 2 (l - {uj/k) 2 fl 2 
[l+p 2 sin 2 a + fi2 g 2)i/2' 



(4) 



In the system defined by F, variables (p, q) are slow and variables (P, 4>) are fast. Slow 
evolution of (p, q) is determined by a system with Hamiltonian ^/eh. This system defines 
a flow on the resonant surface, called resonant flow. The (P, <ft) variables evolve according 
to the subsystem with Hamiltonian Fq. If < b < 1, there is a separatrix surrounding the 
oscillation region on the phase portrait of this subsystem (see Figure ||). If b < or b > 1, 
there is no oscillation region. 

The area of the oscillation region S is a function of the slow variables: S = S(p,q). 
If S(p, q) enlarges along the resonant flow, additional area appears inside the oscillation 



region. Hence, phase points cross the separatrix and enter the oscillation region. This is 
a capture into resonance. A captured phase point leaves a vicinity of the curve I =const, 
H,q =const and continues its motion following approximately the resonant curve. Note, that 
phase points with arbitrarily large initial values of P can be captured provided they are 
close enough to the incoming invariant manifold of the saddle point of the "pendulum- like" 
system. This corresponds to the fact that a particle can be trapped in the mode of surfing 
acceleration even in the case that initially it is far from the resonance. 

The area bounded by the trajectory of a captured phase point in the (P, 0)-plane is an 
adiabatic invariant of the "pendulum-like" system. Hence, if S(p, q) contracts along the 
resonant flow, some phase points leave the oscillation region and leave the resonant zone. 
This is an escape from the resonance. If S monotonically grows along the resonant curve, 
none of phase points leave the oscillation region. In this case, captured phase points stay 
captured forever. 

If the resonant curve is a hyperbola (k 3 < u, @j) or a parabola (k 3 = u,Ho < 0, 0), 
a captured phase point may go to infinity. In this motion the energy of the particle H 
(see ([I])) tends to infinity. Therefore, this motion produces unlimited surfing acceleration 
of particles. This acceleration is possible, if S(p, q) grows as p, q — ► oo along the resonant 
curve. Calculations (see [Q) give the following necessary condition of possibility of unlimited 
acceleration: 



This condition was first obtained in |3J for acceleration by a harmonic wave. In the case of 
perpendicular propagation, it is equivalent to the condition of Katsouleas and Dawson ||, 
also mentioned in [I]]. 

Consider a hyperbolic resonant curve under assumption that @ is valid. At q < 0, 
b(p,q) < (see (|)) and S(p,q) =0. At q = 0, function S(p,q) has a singularity and at 
small positive q it is very large. As q grows along the resonant curve, S(p, q) first decreases 
and then, as q — > oo, S(p, q) — > oo. Hence, at a certain q = q m , function S(p, q) has minimum 
S = S m along the resonant curve. Consider a phase point that is initially captured into the 
resonance at small positive value of q. Let the area bounded by its trajectory be So and 
S > S m . Then in the course of motion along the resonant curve this phase point escapes 
from the resonance. If So < S m , the phase point stays captured forever and undergoes 




(5) 
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unlimited acceleration. This explains Fig. 5 in 0. The number of bounces in this figure 
corresponds to the number of oscillations of the phase point inside the oscillatory region of 
the "pendulum- like" system, performed before the phase point escapes from the resonance. 

The author thanks A.I.Neishtadt for useful discussions. The work was partially supported 
by grants RFBR 00-01-00538 and INTAS 00-221. 
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